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Abstract. We study the critical behavior of Ising quantum magnets with broadly distributed random 
couplings (J), such that P(ln J) ~ | In J|~^~°, a > 1, for large | In J| (Levy flight statistics). For sufficiently 
broad distributions, a < ac, the critical behavior is controlled by a line of fixed points, where the critical 
exponents vary with the Levy index, a. In one dimension, with Oc — 2, we obtaind several exact results 
through a mapping to surviving Riemann walks. In two dimensions the varying critical exponents have 
been calculated by a numerical implementation of the Ma-Dasgupta-Hu renormalization group method 
leading to ac ~ 4.5. Thus in the region 2 < a < a^, where the central limit theorem holds for | In J\ the 
broadness of the distribution is relevant for the 2d quantum Ising model. 

PACS. 75.50.Lk Spin glasses and other random magnets - 05.30. Ch Quantum ensemble theory -75.10.Nr 
Spin-glass and other random models - 75.40.Gb Dynamic properties (dynamic susceptibility, spin waves, 
spin diff'usion, dynamic scaling, etc.) 



1 Introduction 



Quenched, i.e. time independent disorder may have a 
strong influence on quantum phase transitions, which take 
place at zero temperature by varying a quantum con- 
trol parameter, S, which measures the strength of quan- 
tum fluctuations |l||. The critical behavior of random quan- 
tum magnets can be conveniently studied within the 
framework of a renormalization group (RG) scheme intro- 
duced bvMa, Dasgupta and Hu|g] and later applied by 
Fisher llQ and others ||,|,§. In this "energy-space" RG 
method the strongest bonds and transverse fields are suc- 
cessively decimated out and other couplings/fields are re- 
placed by weaker ones generated by a perturbation cal- 
culation. During the RG procedure one keeps track of 
the energy scale, i7, which is the actual value of the 
strongest coupling/transverse field; the average length 
scale, £,, which is the correlation length associated with 
the average equal-time correlation function; and the size 
of the average local moment of a cluster, fj,, which is 
formed after decimating out the couplings between the 
spins. For many systems - during renormalization - the 
disorder grows without limits, thus as J7 — > typically a 
coupling/transverse field is infinitely stronger or weaker 
than the neighboring one and the RG procedure becomes 
asymptotically exact [H. This type of critical behavior is 



controlled by an infinite randomness fixed point (IRFP), 
At an IRFP the physical properties of the system are de- 
termined by the so called rare events, which are samples, 
which occur with vanishing probability, but dominate the 
average behavior. 

At an IRFP the scaling properties are unusual: they 
are fully characterized by three critical exponents, denoted 
by -0, (/) and v and defined as follows P|. At the IRFP 
the relevant time-scale, tr ^ 1/^, and length-scale, £,, are 
related through: 

lni,^e^ (1) 

thus scaling is strongly anisotropic and, as a consequence, 
the average critical dynamical correlations are logarithmi- 
cally slow|jlQ|. The dependence of the average magnetic 
moment per cluster on the change of the energy scale as 
Qq —f f2 is asymptotically given by 



M ^ Mfio/O)]" 



(2) 



from which the anomalous dimension of the bulk magne- 
tization, Xm, in dimension, d, follows as: 



Xm = d — (pip 



(3) 



We note that for bulk spins the average, equal-time cor- 
relations behave as ~ r~'^^"^ at criticality, whereas end- 
to-end correlations involve the corresponding surface ex- 
ponent, x^, which is generally different form Xm- Finally, 
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Table 1. Critical exponents of the RTIM in one and two di- 
mensions. In Id exact results for normal disorder are in the 
first row, whereas in the second row the corresponding expo- 
nents are given for the Levy distribution (1 < a < 2). In 2d 
the differences in the numerical estimates are due to different 
disorder distributions and the truncation procedure used in the 
numerical KG method. 
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the effect of inhomogeneous or correlated disorder have 
aheady been studied in Ref . jlj] and Ref . ||l^ , respectively. 
In both cases the critical properties for strong enough 
perturbations are modified: they are governed by a line 
of IRFP-s, such that the critical exponents are continu- 
ously varying functions of the inhomogeneity/correlation 
parameters. 

In the present paper we are going to release the third 
type of restriction and study the effect of broad disorder 
distributions on the critical properties of random quantum 
magnets. We consider the RTIM with the Hamiltonian 
in Eq.(|^) and keep in mind that at the IRFP it is the 
logarithm of the couplings and the transverse fields which 
follows a smooth probability distribution [pj. Therefore we 
use a parametrization 



the third characteristic exponent of the IRFP, v^ is con- 
nected to the asymptotic behavior of the average correla- 
tion length close to the critical point: 



i-\&\- 



(4) 



Prototype of a disordered quantum system with an IRFP 
is the random transverse-field Ising model (RTIM) defined 
by the Hamiltonian: 



H 



E-^. 



Y.^^< 



(5) 



<ij> 



Here crf , erf are Pauli matrices at site z, the nearest neigh- 
bor coupling constants, J^, and the transverse fields, hi, 
are independent random variables. 

In one dimension (Id), where the topology of the 
lattice is invariant under renormalization, Fisher|3| has 
solved analytically the fixed-point RG equations from 
which the position of the random critical point, 



[ln/i],,-[lnJ], 



0, 



(6) 



and the values of the critical exponents follow. They are 
listed in Table 1. 

In Eq.(||) and in the following [. . .]av is used to de- 
note averaging over quenched disorder. The critical expo- 
nents in Table 1, except ^, can also be derived by an exact 
mapping to a random walk (RW) problem ||ll|. A series of 
numerical studies [|l2| have confirmed the validity of the 
above RG and RW results. 

Also in two dimensions the critical behavior of the 
RTIM is found to be controlled by an IRFPg and the 
corresponding critical exponents, as presented in Table 1, 
have been determined by implementing numerically the 
RG procedure |5|P,0. These RG results are consistent with 
values obtained by quantum Monte Carlo simulations p3| . 

One important question concerning the critical behav- 
ior of random quantum systems is the domain of attraction 
of the IRFP for different type of disorder distributions. 
Using an analogy with the RW problem it is generally 
expected that the critical exponents of random quantum 
systems are universal, provided the disorder is i) spatially 
homogeneous, ii) not correlated, and in) not broadly dis- 
tributed. Problems related to the first two conditions, i.e. 



hi = hQ . 



(7) 



where the exponents, Oij are independent random vari- 
ables. They are taken from a broad distribution, 7r(0), 
such that for large arguments they decrease as, 7r(6') ~ 
16*1^^^". We consider the region a > 1, where the K-th 
moment of the distribution exists for k < a. This type 
of distribution, which comes for in different domains of 
physics and science Ijlg] is usually called Levy flight or Rie- 
mann walk in the discrete version its parameter, a, is 
the Levy index. Throughout this paper we use the fol- 
lowing distribution: 7r(6') — pa{l + 6*)^^^" for 6* > 
and. 7r(6)) = qa{l + |0|)-i-" ior O < 0, p + q ^ 1. In 
one dimension the quantum control parameter in Eq.(|g) 
is given hy S — Inho — {p — q)\nA/{a — l). We also use the 
discretized version (Riemann walk) of the above distribu- 
tion, where 0i = ±1,±2,... etc. and the normalization 
a is replaced by (C(l -I- a) — 1), where ({x) denotes the 
Weierstrass zeta function. 

The structure of the paper is the following. In Section 2 
the Id version of the model is considered and studied ana- 
lytically (through a mapping to surviving Riemann walks) 
and numerically. In Section 3 the 2d model is studied by 
a numerical implementation of the RG procedure. Our re- 
sults are discussed in the final Section, whereas some de- 
tails about persistence properties of Levy flights are pre- 
sented in the Appendix. 



2 The one-dimensional problem 



2.1 Exact results through a mapping to random walks 



The critical properties of the one-dimensional RTIM can 
be conveniently studied through a mapping to a RW 
problem fl^ and a similar procedure works also for the 
random XX- and XY- models p7|. The method, which 
has also been used for inhomogeneous |Q and correlated 
disorder [|l^, is based on an exact expression of the surface 
magnetization of the transverse-field Ising model with L 
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sites @[T): 



TOs(L) 



L-1 I 



1=1 3 = 1 



J, 



-1/2 



(8) 



where the last spin of the chain at Z = L is fixed to the 
state |(t£) = | t)- Before analysing Eq.(||) we cite another 
simple relation |]l^,0, in which the lowest excitation en- 
ergy, e{L), in a finite system with open boundary condi- 
tions is asymptotically related to the surface magnetiza- 
tion in Eq.(p|): 



L~l 



e(L) ~ ms{L)ms{L)hL W -j , 



•h 



(9) 



provided t(L) vanishes faster than L~^. Here Tns{L) de- 
notes the finite-size surface magnetization at the other end 
of the chain and follows from the substitution hj/Jj <-^- 
Hl-j/Jl-] in Eq.(l). 

Now we start to analyse the expression in Eq.(H) and 
look for the possible values of ms{L) using the discrete 
version of the distribution, 71(6*), fix /iq = 1 and taking 
the limit yl — > 00. It is easy to see that for this extreme 
distribution the products in Eq.(||), Hi^i i^j/Jj) 7 take 
three different values: zero, one or infinity, and for a given 
sample nis is zero, whenever any of the products is infi- 
nite, otherwise ms{L) = 0(1). To calculate the average 
surface magnetization one should collect the samples with 
ms{L) = 0(1). Here we use the RW-picture of Ref.[pl| 
and assign to each disorder configuration a random walk, 
which starts at i = at position y — and takes at time 
t = id. step of length 0i with probability ^{Oi). Then for a 
disorder configuration with a finite surface magnetization 
the corresponding RW stays until t = L steps at one side 
of its starting position, 2/{i) > 0, t = 1, 2, . . . , L, in other 
words the RW has surviving character. As a consequence 
[ms(i)]av is proportional to the fraction of surviving RW- 
s, given by the surviving probability, Psurv{t) at t = L. 

For a symmetric distribution, i.e. with p — q — 1/2 
the corresponding RW-s have no drift, whereas for the 
asymmetric case, p y^ q, there is an average bias given by 
Sw = q — p, so that for Sw > 0(< 0) the walk is drifted 
towards (off) the adsorbing wall at y = 0. The bias of the 
RW is proportional to the control-parameter of the RTIM, 
5, as defined in Eq.ra), thus the correspondence between 
RTIM and RW can be generally formulated as: 



[msiS,L)l 



,{Sw,t)\ 



t=L ' 



(10) 



Consequently from the persistence properties of Levy 
fiights, which are summarized in the Appendix, one can 
deduce the singular behavior of the average surface mag- 
netization of the RTIM. 

We start with the finite-size behavior at the critical 
point, (5 = 0, which is given with the correspondences in 
Eqs.(|l^ and ( |AT^ ) as 



K(o,i)], 



1/2 



(11) 



Thus the anomalous dimension of the average surface mag- 
netization, x^, does not depend on the Levy index, a, its 
value is the same as for the normal distribution in Table 
1. 

In the paramagnetic phase, (5 > 0, the corresponding 
Riemann walk has an average drift towards the adsorbing 
site. Consequently its surviving probability in Eq.(A.14) 
and thus the related average surface magnetization of the 
RTIM has an exponentially decreasing behavior as a func- 
tion of the scaling variable, 5L^~^/°' , which is analogous to 
that in Eq.(A.4). Consequently the characteristic length- 
scale in the problem, the average correlation length, ^, 
and the quantum control parameter, (5, close to the crit- 
ical point are related as in Eq.(|4[), however with an a- 
dependent exponent: 



v{a) 



1 



(12) 



Note that v{a) is divergent as a — ^ 1+, which is a con- 
sequence of the fact that the first moment of the Levy 
distribution is also divergent in that limit. In the other 
limiting case, a ^ 2~, we recover the result for the nor- 
mal distribution in Table 1. 

In the ferromagnetic phase, 5 < Q, the correspond- 
ing Riemann walk is drifted off the adsorbing site and, 
as shown in the Appendix, the surviving probability ap- 
proaches a finite value in the large time limit. Conse- 
quently the average surface magnetization of the RTIM 
is also finite in the ferromagnetic phase and for a small \5\ 
it behaves according to Eq.(A.15) as: 



lim 



-M.L)U-W% /3, 



2(a-l) 



(13) 



x^z^, is satisfied. 



Thus the scaling relation, (3^ 

Next, we turn to study the scaling behavior of the 
lowest excitation energy starting with the expression in 
Eq.(y). Here we note that in a given sample the existence 
of a very small gap is accompanied by the presence of local 
order. Thus in a sample with a low-energy excitation one 
can find a strongly coupled domain (SCD) of size I, where 
the coupling distribution follows a surviving walk charac- 
ter. Consequently in Eq.(^ m^ = 0(1) and nis = Oil) 
and one gets the estimate: 



i-i 

n 



J. 



exp 



(-ItrHJ/h)) ■ (14) 



Here Itr measures the size of transverse fiuctuations of the 
corresponding surviving Levy fiight of length t = I and 
ln(J//i) denotes an average value. At the critical point, 
S ^ Sw — 0, the surviving region of the walk and thus the 
SCD in the RTIM exte nds over the volume of the sample, 
I - i, and from Eq.(|0|), kr ^ L^^", so that Eq.^ 
leads to: 



\ne{L) r^ L"^', ^j^l/a 5 = 0; 



(15) 



The critical exponent, ^l' defined in Eq.(nf) is a continuous 
function of the Levy index which takes its value for the 
normal distribution in Table 1 when a ^ 2~ . 
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In the paramagnetic phase, (5 > 0, there are stiU re- 
ahzations with a smaU energy gap, the scahng form of 
which in a finite system of size, L, can be estimated by 
the foUowing reasoning. For a Levy flight the probabil- 
ity of a large transverse fluctuation, Itr, is given from 
Eq.(A.8) as p{ltr,L) ~ Liltr) " , thus its character- 
istic size can be estimated as Itr ~ L^/'^^+") from the 
condition p{ltr,L) = 0(1). Consequently the lowest gap 
has the scaling form: 



3.0 IT 



lne(i)-ii/(i+"), 6>0, 



(16) 



which implies a logarithmically broad gap distribution 
even in the Griffiths phase. This is in contrast to the be- 
havior with the normal distribution, where Itr ~ lnLpl|, 
and the scaling form of the gap is in a power-law form, 
e(L) ~ L~^, where the dynamical exponent, z, is a con- 
tinuous function of the quantum control parameter. 

In the remaining part of this subsection we discuss the 
probability distribution of the surface magnetization. Let 
us remind that, at the critical point, the average surface 
magnetization is determined by the so-called rare events, 
which are samples having rus — 0{1). The typical sam- 
ples, however, which are represented by non-surviving ran- 
dom walks, have a vanishing surface order in the thermo- 
dynamic limit. For a large but finite system of size L, 
ms{L) is dominated by the largest product, Ylj{hj/Jj), 
in Eq.dsl) so that 



liirrisiL) ^ e{L) ^ —ltrln{h/J) 



(17) 



where e{L) is the value of the gap in the dual system, 
i.e. where the fields and the couplings are interchanged, 
/ij <-> Jj. Since at the critical point the system is self-dual 
we have from Eq.(n3) lnms{L) ^ L^ and the appropriate 
scaling variable is 



1 flnm. 



(18) 



Assuming that the scaling function behaves as p{y) ~ 
\y\°' for small |y|, then the average magnetization is given 
by [m,]^^ = l/L'^/dm,p(lnm,/L'^) ~ £"(1+")'^, which 
implies a — a/2 — 1, to recover the exponents in Eqs.([ll|) 
and (ITsI). These relations have been checked through a 
numerical evaluation of Eq. (^ for 50000 samples of large 
systems with size up to L = 2048. As seen in Fig.l there 
is a very good agreement between the scaling and numer- 
ical results. In the paramagnetic phase, 5 > 0, close to 
the critical point according to Eq.^) the typical surface 
magnetization behaves asymptotically as 



InTOs(L) 



-5L^ 



-i/e 



typ 1 



(19) 



where the typical correlation length, S,typ, diverges at the 
critical point as: 



6. 



i-^r 



l^typ — 1 



(20) 



Note that the characteristic exponent, vtypAs independent 
of a, but still satisfies the scaling relationg, vtyp = ^^(1 — 




Iln m l/L ■ 



Fig. 1. Rescaled distribution of the logarithmic surface magne- 
tization for L = 512(A), 1024(o) and 2048(o) at a = 1.5. The 
full line represents the scaling expression for small arguments, 
where the only fitting parameter is the amplitude. 



2.2 Bulk magnetization: a numerical renormalization 
group study 

The critical behavior of the bulk magnetization of the Id 
RTIM is not related to the properties of a homogeneous 
RW, but it can be calculated from the so called average 
persistence properties of a Sinai walk, i.e. a random walk 
in a random environment [pO[ . This procedure has already 
been used to calculate the magnetization scaling dimen- 
sion, Xm, of the RTIM with correlated disorder ||l5|. Here 
we use another numerical method to study Xm for Levy- 
type disorder, which method is based on a numerical im- 
plementation of the Ma-Dasgupta-Hu RG procedure. The 
RG-equationsH and their numerical use are well docu- 
mented in the literature yJ6[0]. Here we use the finite-size 
version of the method||], in which we start with a large fi- 
nite ring of size L with random couplings and perform 
the decimation procedure until the last spin. The last 
log-coupling (transverse field) sets the log-energy scale, 
r — — log(J7/i7o)j which at the critical point scales from 
Eq.(l) as F ~ L'^ , whereas the critical cluster moment in 
Eq.(2), associated with the last remaining cluster, scales 
as /i ~ L*^"^. Repeating the calculation for several realiza- 
tions of the disorder the critical exponents can be deduced 
from the appropriate scaling functions. 

We start with an analysis of the log-energy distribu- 
tion, which should be a function of the scaling variable 
r/L"^ for different finite systems. As shown in Fig. 2 there 
is an excellent data collapse using our analytical result 
t/i = 1/a as given in Eq.([15|). Next, we analyze the dis- 
tribution of the cluster moment /i, which is a function of 
the scaling variable fi/L'^^ , from which the magnetization 
scaling dimension, Xm, follows through Eq.dS). From the 
optimal data collapse, as presented in Fig. 3, we obtain 
Xm = 0.22 for a = 1.5, which is definitely larger than that 
of the system with normal disorder, as given in Table 1. 
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Fig. 2. Rescaled distribution of the logarithmic energy scale, 
r, for different finite lengths, L = 32(o),64(+) and 128(*) 
at a = 1.5. The best collapse of the data is obtained by the 
analytical result tp = 1/a. 



Repeating the calculation for the average magnetization, 
771 ((5) = [/i/L]av, the appropriate scahng variables outside 
the critical point are mL^™ and LS'^. As shown in Fig. 4 in 
terms of the scaled variables we obtain a very good data 
collapse, using Xm from Fig. 3 and the analytical expres- 
sion in Eq.(n3) for v. 



1000 



J 
H 

oT 




Fig. 3. Rescaled distribution of the cluster moment for dif- 
ferent finite lengths, L = 32(-|-), 64(o) and 128(o) at a = 1.5. 
The optimal collapse of the data gives Xm = 0.22. 



For a systematic study of the bulk magnetization scal- 
ing index we have determined the average critical mag- 
netization for different finite systems and then calculated 
an effective exponent, Xm{L), by two-point fit comparing 
m{L) and m{2L). As shown in Fig. 5 the effective expo- 



Fig. 4. Scaling plot of the average bulk magnetization for 
diff'erent finite systems at a = 1.5. 



nents have a weak size dependence, so that one can have an 
accurate extrapolated value for L —>■ oo, which clearly de- 
pends on the value of the Levy index, 1 < a < 2. These ex- 
trapolated magnetization exponents are presented in Fig. 
6, where for a > 2 Xm is expected to be a independent, 
however the corrections to scaling are strong, especially 
around the cross-over value ac = 2. At this point we note 
that in parallel to the RTIM we have also calculated the 
magnetization scaling dimension of the random quantum 
Potts chain for q = 3. 



025 



0.24 



022 



021 



0.2 



0.19 



0=1.0 




O05 0.1 

1/L 



Fig. 5. Effective magnetic exponents for a — 1.5 and a — 1. 
The full lines connecting the data for the RTIM are guide to 
the eye. For a comparison the same quantities for the q = 3- 
state random quantum Potts chain are also presented, where 
the data points are connected by dashed lines. 



As can be seen in Fig. 5 the extrapolated exponents of 
the Ising and Potts chains are very probably the same for 
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the same value of the Levy index, a. This result completes 
the universality of the two models as obtained before an- 
alytically for normal disorder ||21|. 

We close this section by a study of the average magne- 
tization in the ferromagnetic phase, S < 0, where the aver- 
age cluster moment scales with the size of the system, L. 
As a consequence the average magnetization approaches 
a finite limit as L — > oo, which close to the critical point 
behaves asymptotically as m{S) ^ \S\^. As shown in Fig. 7 
the results for different finite systems converge to a power- 
law form where the critical exponent, /3, satisfies the scal- 
ing relation, /? = Xmi^. 
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Fig. 6. Bulk magnetization exponent for different values of 
a. The full line represents Xm for normal disorder, the dashed 
line is guide to the eye. 



3 The two-dimensional problem 

In 2d the critical properties of the RTIM with normal dis- 
order is controlled by an IRFP, so that the Ma-Dasgupta- 
Hu RG procedure provides asymptotically exact critical 
properties. We expect that these properies of the RTIM 
stay valid for broad distributions, too, and we study the 2d 
problem using the RG method. There are, however, sev- 
eral limitations for the numerical implementation of the 
method. First, as in Id, the systems under consideration 
have a finite spatial extent and the number of disorder 
realizations which can be studied is also finite. In 2d an 
extra problem is caused by the fact that the topology of 
the lattice is not invariant under the RG transformation 
and as a result new couplings between remote spins are 
generated. 

To avoid these type of difficulties a truncation proce- 
dure has been introduced in RefJ5|, such that a class of 
generated interactions, which are expected to cause very 
small errors, are neglected. In Ref. [||]], where the finite- 
size version of the RG method is used, all the generated in- 
teractions were kept for moderately large finite systems. In 



Fig. 7. Average bulk magnetization in the ordered phase as 
a function of the quantum control parameter in a log-log plot 
(a = 1.5). For small S the finite size results converge to the 
dashed line with a slope oi fi — 0.66. 



Ref. [0] a selection condition in the finite-size RG method 
has been introduced by showing that many of the gen- 
erated new couplings are "dead" in the sense that they 
are not decimated out in later steps, so that they can be 
omitted without causing any error in the renormalization. 
With this observation one could consider larger finite sys- 
tems and at the same time the necessary computational 
time has considerably reduced. 

In the present paper we apply the finite-size RG 
method supplemented by the selection condition. In this 
way we could treat systems on the square lattice with lin- 
ear size up to L = 128 and we considered typically 10000 
realizations and some 1000 samples for the largest sys- 
tems. 

In 2d the position of the critical point of the system is 
not known by self-duality, therefore we used the following 
numerical procedure for its determination. First we cal- 
culated the average magnetization, m{L,ho), in a finite 
systems of size L and then defined the scaling function, 
5l(/io), as the ratio: 



m{L,ho) 
m(L/2,/io) 



9L{ho) 



(21) 



In the ferromagnetic phase, /iq < /ic, in the thermody- 
namic limit m(L, /ig) does not depend on L, consequently 
lim^^oo ffL(^o) = 1- On the other hand in the paramag- 
netic phase, /iQ > he, using the exaniple of the surface 
magnetization of the Id RTIM in Eq.(@), m{L,hQ) is ex- 
ponentially small in L, thus the scaling function has a 
vanishing limiting value. In between at the critical point, 
ho = he, where m{L, he) '^ L~^"' we have a finite limiting 
value liuiL^oo gL^hc) — 2^^™. Consequently calculating 
51, (/iq) for a series of sizes the position of the limiting 
crossing points defines the critical point of the system, 
whereas the abscissas of the crossing points are related to 
the magnetization scaling dimension, Xm- This procedure 
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is illustrated in Fig. 8 on the example of the Id model with 
a = 1.5, where the previously calculated critical proper- 
ties of the system are accurately reobtained. 



* * L=2048 

D B.L=1024 

L=512 
L=256 




Fig. 8. Finite-size estimates of the critical point and the mag 
netization scaling dimension of the Id model with a = 1.5. 



value can be estimated as ac ~ 4.5. The magnetization 
scaling dimension, x™, or more precisely (l)ip — 2 — Xm, 
is also related to the finite-size behavior of the average 
cluster moment at the critical point: [/x^Jav ~ L'^'^ , which 
is illustrated in Fig. 11 for a = 3. From the slope of the 
curve in a log-log plot one obtains an estimate of (j)ip = .48, 
which is in agreement with the results shown in Fig. 10. 
Similar agreement is found for other values of a, too. 
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Fig. 9. The same as in Fig. 8 for the 2d model with a = 3. 



In 2d a similar plot of the scaling functions for a = 3 
is shown in Fig. 9, where the crossing point is obtained 
with a reasonably small error. Repeating this calculation 
for various a-s we have obtained estimates for the critical 
points and the magnetization exponents, which are pre- 
sented in Fig. 10. As seen in this figure Xm is monotonously 
decreasing with a < ac, whereas for a > ac Xm stays ap- 
proximately constant and this value corresponds, within 
the error of estimates, to that of the 2d RTIM with normal 
disorder, as given in Table 1. ^From Fig. 10 the cross-over 



Fig. 10. Magnetization scaling dimensions in 2d for different 
values of the Levy-index. In the inset the estimated critical 
points are presented. The line connecting the data points is a 
guide to the eye. 
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Fig. 11. Average critical cluster moment as a function of size 
in a log- log plot {a — 3). The slope of the line corresponds to 
(pip = 0.48. 

Our final investigation concerns the scaling behavior of 
the distribution function of the log-energy-scale, Pl(/^), 
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at the critical point. As in the Id problem the critical 
exponent, ip, can be obtained from an optimal data col- 
lapse in terms of the scaling variable F/ L^ . In Fig. 12 we 
present for a = 3 the rescaled accumulated probability 
distribution function, where a satisfactory data collapse 
is obtained for i/; — 0.8. Similar estimates for other a- 
s are collected in Fig. 13. As for the Xm exponent ip is 
a monotonously decreasing function for a < ac '^ 4.5, 
whereas for a > etc it is approximately constant and this 
value corresponds to that of the model with normal dis- 
order as given in Table 1. 




r/L 



Fig. 12. Rescaled accumulated probability distribution func- 
tion of the log-energy-scale at a = 3. Satisfactory data collapse 
is obtained with the exponent ip — 0.80. 




Fig. 13. The critical exponent i/> as a function of the Levy 
index for the 2d model. The line is a guide to the eye. 



4 Summary 



In this paper the effect of a broad disorder distribution 
on the critical behavior of the RTIM is studied in one 
and two space dimensions. The broadness of the disorder 
distribution becomes relevant, if the Levy index is lowered 
below a critical value, Uc- In the region of 1 < a < etc the 
critical exponents of the IRFP-s are continuous functions 
of a and for a > ac they are the same as in the model 
with normal disorder. 

In Id we obtained ac — 2, in close analogy with ran- 
dom walks, where the central limit theorem is valid for 
a > 2. This analogy is more than a simple coincidence, 
since the Id RW and the quantum Ising spin chain are 
related through an exact mapping |22[| , which then has the 
same requirement for the relevance- irrelevance conditions. 
In 2d, where such type of mapping does not exist, the lim- 
iting value is found to be approximately around ac ~ 4.5, 
thus in the region of 2 < a < ac the broadness of disorder 
is relevant for the RTIM, whereas it is irrelevant for the 
RW. 

Outside the critical point, in the so-called Griffiths- 
phase ||23[ some physical quantities of random quantum 
Ising magnets (linear and non-linear susceptibility, auto- 
correlations, etc.) are still singular and in the presence 
of normal disorder these Grifhths-McCoy singularities are 
characterised by the dynamical exponent, z{6), which is 
a function of the distance of the critical point. For broad 
distribution of the disorder the probability of rare events is 
enhanced, consequently the GrifRths-McCoy singularities 
become stronger. In Id we have shown by scaling argu- 
ments leading to Eq.(|l6|) that the dynamical exponent is 
formally infinite in the whole Griffiths-phase. 

Considering other magnetic models we expect that in 
Id the critical properties of the random XX, XXX and 
q-state quantum Potts models are controlled by the same 
IRFP as that of the RTIM and this is parametrized by 
the Levy index, a. This result for the XX- model is a con- 
sequence of a known mapping, [|iyi 7]] whereas for the two 
other models the RG decimation equations are asymptot- 
ically equivalent to that of the RTIM|^,|l|. This equiva- 
lence has been explicitly tested here for the Potts model 
with q = 3. On the other hand in two dimensions such an 
equivalence can be expected only for the random quantum 
Potts model, since the random XX- and XXX-models have 
no IRFP-s in higher dimensions fiO . 
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Appendix: Surviving probability of Levy flights probability for Sn > vn, irrespectively from the previous 

steps, is denoted by P{n,v) and the latter is given by: 
Consider the following sum 

P{n,v)^ / p{S,n)dS . (A.9) 

j=i Between the generating functions: 



where the independent random variables, Xj, follow the 
same broad probability distribution, 7r(x), which asymp- 
totically behaves as 



nix) ~,^oo p.T-(i+") 



tt{x) 



q\x\ 



~(l+a) 



P{z,v 



n>0 






(A.IO) 



(A.2) 

with 1 < a < 2 and we are interested in the probability 
distribution of 5„, p{S, n), for large n. According to exact 
results pJ there exists a limit distribution, p{u)du, in term 
of the variable, u = S'„/Z„ — c„, as n —> oo. Here the 
normalization: 



Ir. 



1/a 



(A.3) 



is the transverse fluctuation of the walk, if we interpret 
n — t as the (discrete) time and Sn=t as the position of 
the walker in the transverse direction. The second normal- 
ization is given by 



there is a useful relation due to Sparre Andersenpsf: 

F{z,v) = exp[P(z,w)] , (A.ll) 

which has been used recently in Ref||26|]. 

In the zero velocity case, v — 0, which is equivalent to 
the symmetric walk with S^ = 0, we have P(n,0)=l/2. 
Consequently P{z,0) = -iln(l - z) and F(z,0) = (1 - 

z)~^/^, from which one obtains for the final asymptotic 
result: 



= -ni-i/"<5„ 



(A.4) 



Psurv{t,0)= F{n,0)\ 



r 



1/2 . (A.12) 



where with 6^, = — (x) we define the bias of the walk . For 
a small Sw one gets from the combination in Eq. ( A.4 ) the 
scaling relation between time and bias as: 



t^\S^r^"\ v{a) 



a — 1 



(A.5) 



For a symmetric distribution, when tt^x) is an even func- 
tion, thus p = q and S^ = 0, we have for the limit distri- 
bution: 



Note that the persistence exponent, 9 ~ 1/2, is indepen- 
dent of the form of a symmetric probability distribution, 
7r(x), thus it does not depend on the Levy index, a. 

For V > 0, i.e. when the allowed region of the particle 
shrinks in time the correction to P{n, 0) = 1/2 has the 
functional form ^ P( n, v) = 1/ 2 — q (c), with c — vv}~^l°'. 
Evaluating Eq.(A.9) with Eq.( |A.7D one gets in leading or- 
der, P{n,v) = y= g^(a) -F 0(P), with A{a) = r(l -f 
l/a)/7r. Then P(z, v) - P{z, 0) ~ A{a)v X;„>i z'^n,-^/"' is 
singular around z -^ 1~ as^^ (1 — z)^(^^^/"\ consequently 



1 r°° 

p{u) = L^^oiu) = — J e"^-''-l'=l°dfc, (A.6) F{z,v) ^ {1 - z)-^/' 



exp 



which has an expansion around u ^ 



u- „ /2fc + 1 



1 oc> 2k 



(A.7) 



(A.13) 
in leading order and close to z = 1 . Here the second fac- 
tor gives the more singular contribution to the surviving 
probability, which is in an exponential fornr: 



P.. 



,it,S^) - P(n,w)|„^5^,„^t 



and for large u it is asymptotically given by 

L„ o(u) = 1m-(i+")p(1 + a) sin(7ra/2) , (A.8) 

TT 

where r{x) denotes the gamma function. 

Consider next the surviving probability, Psurv{i^5w)i 
which is given by the fraction of those walks, which have 
not crossed the starting position until t — n, thus Si > 
for z = 1, 2, . . . , n. For a biased walk, with < |5^| ^ 1, 
the asymptotic behavior of Psurvin^^w) is equivalent to 
that of a symmetric walk {5w — 0) but with a moving 
adsorbing boundary site, which has a constant velocity of 
V — 5w For this event, with Si > vi for i = 1,2, ... ,n, the 
surviving probability is denoted by F{n,v), whereas the 



t 



-1/2 



exp 



-const(5,„i^ ^/"^ 



,(A.14) 



For V < 0, i.e. when the allowed region of the par- 
ticle increases in t ime we co nside r the large \v\ limit 
and write Eq. ( 
B{a)£ 



with Eq.(A 



-0{c- 



as P{n,v) ~ 1 — 
with B{a) = r{l + a) sm{TTa/2)/TTa. 
Then, in the large |v| limit P{z,v) — — ln(l — z) — 
B{a)\v\~" '^^y^ z"n~°', where the second term is con- 
vergent even at z = 1. As a consequence the surviving 
probability remains finite as n — > cx) and we have the re- 
sult, F{n,v) ~ 1 — const|u|~", for \v\ ^ 1. For a small 
velocity, < |u| <^ 1, we can estimate F{n,v) by the 
following reasoning. After n — Uc steps the distance of 
the adsorbing site from the starting point, ys = vric, will 
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exc eed the size of transverse fluctuations of the walk in 
Eq(|A^), kr ~ n^", with ric ^ \v\-'''-"\ Then the walker 
which has survived until nc-steps with a probability of 

— 1/2 

Uc , will survive in the following steps with probability 
0(1). Consequently 

lim Psurvin,S^) ^ lim F(n,v)|,„^^ „^, ~ I^^T^")/' . 

(A.15) 
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